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(57) ABSTRACT 

A technique for landing aircraft using an aircraft landing 
accident avoidance device is disclosed. The technique 
includes determining at least two probability distribution 
functions; determining a safe lower limit on a separation 
between a lead aircraft and a trail aircraft on a glide slope to 
the runway; determining a maximum sustainable safe 
attempt-to-land rate on the runway based on the safe lower 
limit and the probability distribution functions; directing the 
trail aircraft to enter the glide slope with a target separation 
from the lead aircraft corresponding to the maximum sustain- 
able safe attempt-to-land rate; while the trail aircraft is in the 
glide slope, determining an actual separation between the 
lead aircraft and the trail aircraft; and directing the trail air- 
craft to execute a go-around maneuver if the actual separation 
approaches the safe lower limit. Probability distribution func- 
tions include runway occupancy time, and landing time inter- 
val and/or inter-arrival distance. 

20 Claims, 13 Drawing Sheets 
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METHOD AND DEVICE FOR LANDING 

AIRCRAFT DEPENDENT ON RUNWAY 
OCCUPANCY TIME 

CROSS REFERENCE TO RELATED 5 

APPLICATIONS 

This application claims the benefit of U.S. Provisional 
Patent Application No. 61/077,230, filed Jul. 1, 2008, which 
is herein incorporated by reference. 10 

STATEMENT REGARDING FEDERALLY 
FUNDED RESEARCH 

This invention was made in part with government support 15 
under Grant Numbers NASA Grant 200932, 200833, Nasl- 
02117, and NSF 200621 awarded by the National Aeronau- 
tics and Space Administration. The Government may have 
certain rights in the invention. 

20 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 shows an example of a flow diagram of a method for 
landing aircraft. 

FIG. 2 illustrates an exemplified pictorial diagram of an 25 
example approach and landing system. 

FIG. 3 shows another example of a flow diagram with 
additional aspects of determining the maximum sustainable 
attempt-to -land-rate. 

FIG. 4 shows an example of a block diagram of an aircraft 30 
landing accident avoidance device. 

FIG. 5 shows exemplified instructions for landing an air- 
craft, where the instructions are stored in a computer readable 
storage medium. 

FIG. 6 shows an example of a bird’ s eye view over the glide 35 
slope and go -around procedures on the glide slope and the 
rates. 

FIG. 7 shows an example of the distributions of ROT and 
LTI. 

FIG. 8 shows an example of total P(GA), P(LTI<ROT), and 40 
dp/doo for WV safe threshold of 65 s. 

FIG. 9 shows landing/h (left) and P{GA} (right) vs. 
attempts/h (oo). 

FIG. 10 shows ^(oo) and g(oo;r) for r=l , 2, and 4 from top to 
bottom, respectively. 45 

FIG. 11 shows ES(co;R,C) of one hour peak period landing 
operations. 

FIG. 12 shows g(oo;r) for different values of cost to benefit 
ratio, where optimal throughput increases as r decreases. 

FIG. 13 shows optimal throughput in r for safe WV thresh- 50 
olds of 65, 70, and 75 s. 

DETAILED DESCRIPTION OF THE 
EMBODIMENTS 

55 

The following notations for input parameters may be used 
throughout the present application. 

LTI^ +1 refers to landing time interval between aircraft k 
and k+1 measured at the runway threshold in seconds and 
assumed to have the lower limit L. 60 

ROT^ refers to runway occupancy time of aircraft k mea- 
sured in seconds. 

R refers to dollar benefit of one successful landing. 

C refers to expected average cost of a go-around or unsuc- 
cessful landing. 65 

x 0 refers to minimum wake vortex (WV) safe separation of 
successive aircraft given in seconds. 


2 

DPI refers to decision point 1 ; nautical mile (nmi) distance 
from threshold where pilot/controller decides whether to 
execute go-around procedure to avoid simultaneous runway 
occupancy (SRO). This is officially referred to as decision 
height. 

DP2 refers to decision point 2; nmi distance from threshold 
where pilot/controller decides whether to execute go-around 
procedure to avoid hazardous wake vortex encounter. 

The following notations for decision variables may be used 
throughout the present application. 

co refers to landing attempts per hour, i.e., flow rate through 
the glide slope, and oo=3 600/mean (LTI) 

A, refers to arrival rate to TRACON or, equivalently, the 
runway throughput rate, landing per hour. 

p refers to probability of go-around P{GA}. 

Reference will now be made in detail to several embodi- 
ments, examples of which are illustrated in the accompanying 
drawings. The embodiments relate generally to a method and 
system for landing aircraft on a runway, and more particularly 
to determining if there is too large a risk of wake vortex or 
simultaneous runway occupancy, and executing a go-around 
maneuver, and otherwise landing the aircraft. Wherever con- 
venient, the same reference numbers will be used throughout 
the drawings to refer to the same or like parts. Although some 
of the following embodiments are described in the context of 
a controller on an aircraft, the disclosed embodiments are not 
limited to being used as a controller on an aircraft and may be 
applied to other types of configurations, including, for 
example, on a distributed system between two aircraft, at a 
control tower for a runway, etc. 

FIG. 1 shows an example of a method 100 for landing 
aircraft on a runway 260. Complimenting FIG. 1, FIG. 2 
illustrates an exemplified pictorial diagram of an example 
approach and landing system 200. To perform the method, an 
aircraft landing accident avoidance device 400 may be used. 
Generally, the method comprises determining at least two 
probability distribution functions 102, where one of the prob- 
ability distribution functions (PDFs) is runway occupancy 
time (ROT), and where another is (1) landing time interval 
(LTI); (2) inter-arrival distance (IAD); or (3) a combination 
thereof; determining a safe lower limit on a separation 
between a lead aircraft 220 and a trail aircraft (also referred to 
as follow aircraft) 230 on a glide slope (also sometimes 
referred to as glide path) to the runway 260, 105; determining 
a maximum sustainable attempt-to -land-rate on the runway 
260 based on the safe lower limit 110; directing the trail 
aircraft 230 to enter the glide slope with a target separation 
425 from the lead aircraft 220 corresponding to the maximum 
sustainable attempt-to -land-rate 115; while the trail aircraft 
230 is in the glide slope, determining an actual separation 436 
between the lead aircraft 220 and the trail aircraft 230, 120; 
and directing the trail aircraft 230 to execute a go-around 
maneuver 240, 250, 445 if the actual separation 436 
approaches or is determined to be approaching the safe lower 
limit 125, 130. If the actual separation 436 is not approaching 
or determined to be approaching the safe lower limit, then the 
trail aircraft 230 is directed to attempt a landing 135. 

Multilateration systems have been installed in some air- 
ports, such as Detroit Metropolitan Wayne County airport 
(DTW), to provide reasonably accurate time-position esti- 
mates of all transponder-equipped aircraft (a/c) operating in 
the airport vicinity in all weather conditions. These data can 
be used to obtain samples of landing process variables, such 
as LTI between successive aircraft to the runway threshold, 
IAD between two successive aircraft at the moment that the 
lead aircraft crosses the runway threshold, and ROT. ROT is 
the length of time required for an arriving aircraft to proceed 
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from over the runway threshold to a point clear of the runway. 
The present invention considers LTI, IAD, and ROT as (ran- 
dom) variables. 

Lower limit (L) (also sometimes referred to as safe lower 
limit) is the minimum time for an LTI distribution. L may be 
based on historical data. LTI is the time it takes between one 
aircraft landing and another aircraft landing. When multiple 
aircrafts are landing on a same path (e.g., runway, landing 
strip, etc.), there needs to be sufficient separation time 
between two or more aircrafts to allow for safe landings. In 
essence, as one embodiment of the present invention, a safe 
lower limit may be determined to be at about 40 seconds. 
Lead aircraft 220 refers to the aircraft that is landing first. 
Trail aircraft 230 refers to the aircraft that is to land immedi- 
ately after the lead aircraft 220. Both aircraft can be any type 
of aircraft. For example, aircraft may be any jet and/or pro- 
peller aircraft, any commercial aircraft (e.g., Delta Airlines, 
Continental Airlines, American Airlines, British Airways, 
FedEx, etc.), private aircraft (e.g., Marquis Jet, etc.), military 
aircraft (e.g, Bl, C-17, C-130, F-16, F-18, KC-135, etc.) 
Sizes of the aircraft may also vary from small (e.g, Cessna, 
Learjet, etc.) to medium (e.g, MD- 80/88, Boeing 737-300/ 
400/500 series, Airbus 320 series, etc.), to large (Boeing 757, 
Boeing 767, Airbus 310 series, etc.) to heavy (e.g, Boeing 
747, Boeing 777, Airbus A300, Airbus 380, etc.). 

Glide path 210 refers to an approach path (which may be 
designated) to a runway 260 or other landing surface for the 
aircraft. Glide path 210 may include the final descent path for 
aircraft. Aircraft may enter the glide path 210 upon approach 
or from a holding pattern to an airport. In one embodiment, 
glide path 210 may include a starting point in three dimen- 
sional space, and an ending point in three dimensional space, 
such as a point on a runway 260 and the path between the 
starting and ending points. To monitor aircraft into and 
through a glide path 21 0, one or more operators from a ground 
control system 270 (such as a control tower, ground based 
sensors (e.g, radar, etc.), etc.) may be needed. 

Throughout the present invention, runway 260 is defined as 
a single runway or at least two parallel runways that are at 
least about 2500 feet apart from each other. 

PDFs and the lower limit may be determined in any order. 
For instance, PDFs may be determined first, followed by the 
lower limit. Alternatively, the lower limit may be determined 
first, followed by the lower limit. As another example, both 
may be determined simultaneously. 

Once the PDFs and the lower limit have been determined, 
a maximum sustainable attempt-to -land-rate can then be 
determined 110. Based on the safe lower limit, the maximum 
sustainable attempt-to -land-rate is the target rate at which 
multiple aircraft is allowed to safely land on the same path- 
way (i.e, runway, etc.). The units for such target are time- 
based, such as be per quarter hour, per half hour, per hour, etc. 

As indicated in FIG. 3, the maximum sustainable attempt- 
to -land-rate includes determining a maximum sustainable 
safe throughput 305, a maximum sustainable net economic 
output 310, or a combination of the two. Sustainable means 
average. As described below, the maximum sustainable safe 
throughput 305 is determined using a probability of executing 
the go-around maneuver 240, 250, 445 and the maximum 
sustainable attempt-to-land-rate. In addition, as described 
below, the maximum sustainable net economic output 310 
(e.g, optimality, profit, benefit, gain, surplus, etc.) is deter- 
mined using a ratio premised on the cost-benefit ratio of 
go-around maneuvers 240, 250, 445 to landings. 

Go-around maneuvers 240, 250, 445 may be any maneuver 
by trail aircraft 230 to fly out of the glide path 210. There may 
be a multitude of reasons why a trail aircraft 230 may execute 
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a go-around maneuver 240, 250, 445. For instance, there may 
be the wake vortex created by the lead aircraft 220 may be so 
large and/or long that it endangers the safety of the trail 
aircraft 230. Another reason may be the trail aircraft 230 may 
5 be entering the glide path 2 1 0 at too fast of a speed that it will 
not have sufficient space to land safely and avoid the lead 
aircraft 220. A further reason may be that one or a flock of 
birds, which may get sucked into the engine(s), may suddenly 
enter the glide path 210. 

10 For whatever reason the trail aircraft 23 0 needs to execute 

a go-around maneuver 240, 250, 445, the go-around maneu- 
ver 240, 250, 445 works in this way: the trail aircraft 230 first 
leaves the glide path 210. After clearing the danger zone and 
15 reentering a safe flight, the trail aircraft 230 may circle 
around. It may also reenter the same holding pattern or enter 
a new holding pattern. Whichever holding pattern the trail 
aircraft 230 has entered, the trail aircraft 230 may then be 
instructed to reenter the glide path 210 and attempt another 
20 landing. 

With the lower limit and the maximum sustainable 
attempt-to-land-rate at hand, the trail aircraft 230 may be 
directed to enter the glide slope with a target separation 425 
from the lead aircraft 220 corresponding to the maximum 
25 sustainable attempt-to-land-rate. While the trail aircraft 230 
is in the glide slope, the actual separation 436 between the 
lead aircraft 220 and the trail aircraft 230 needs to be deter- 
mined. The actual separation 436 is the time that it takes for 
the trail aircraft 230 to reach a specific coordinate that the lead 
30 aircraft 220 once occupied. The specific coordinate can be 
viewed as a marker. If the actual separation 436 approaches 
the safe lower limit, the trail aircraft 230 is directed to execute 
a go-around maneuver 240, 250, 445. “Approaches” means 
the time (e.g, 1 second, 5 seconds, etc.) it takes to reach the 
35 safe lower limit by taking into consideration the pilot’s reac- 
tion time and the response time of the plane from a mechani- 
cal and/or sensor perspective based on the pilot’s reaction. 

Moreover, the timing of when to direct the trail aircraft 230 
to execute a go-around maneuver 240, 250, 445 depends the 
40 reaction time of the trail aircraft 230. Reaction time includes 
the time for the pilot to react, computer response time, time 
for mechanical parts, instruments, and sensors to respond to 
the pilot’s actions, communication, etc. 

As an analogy, the target separation 425 and actual sepa- 
45 ration 436 can be viewed from the perspective of tailgating. 
As a rule of thumb, to determine if one is tailgating on a 
highway, drivers are suggested use a marker on the road (e.g, 
a tree, bridge, sign, etc.) and count 2 seconds from when the 
vehicle in front passes the marker and the vehicular driver 
50 passes that same marker. If the vehicular driver is under 2 
seconds, the vehicular driver is beyond the safe lower limit 
and target separation with respect to braking distance. Thus, 
the actual separation (the distance between the vehicular 
driver and the vehicle in front) is smaller than the target 
55 separation. This shorter distance means that the vehicular 
driver should slow down to increase the distance between the 
vehicular driver and the vehicle in front. 

Along the glide path 210, there may be one or more wake 
vortex decision points where one or more aircrafts may 
60 execute a go-around maneuver 250, 445. As such, the go- 
around maneuver 250, 445 begins at a wake vortex decision 
point. If at any time, such an aircraft misses or bypasses one 
wake vortex decision point, or if such aircrafts fails an 
attempted go-around maneuver 250, 445 the glide path 210 
65 may have additional wake vortex decision points to allow 
execution of go-around maneuver 250, 445 or another 
attempt at a go-around maneuver 240, 445. 
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The wake vortex decision point may be a set of points on 
the glide path 210. For example, the range can be within five 
to twenty nautical miles from the start of runway 260. Each 
wake vortex decision point may be determined as (1) a dis- 
tance measured from the beginning point of the runway 260 to 5 
a point in the glide path 210, (2) a time measured from runway 
260 while being in the glide path 210, (3) a distance and/or 
time in the glide path 210 between the lead aircraft 220 and 
the trail aircraft 230, (4) a distance and/or time in the glide 
path 210 between the trail aircraft 230 and the simultaneous 10 
runway occupancy (SRO) decision point (which is also a 
point or range of points along the glide path 210), or (5) any 
combination of these. Other embodied factors that may be 
considered include the difference in velocity between lead 
aircraft 220 and trail aircraft 230, as well as any expected 15 
change in velocity between lead aircraft 220 and trail aircraft 
230. 

Wake vortex (also generally known as wake turbulence) is 
turbulence, sometimes in the shape of a vortex, that forms 
behind an aircraft as it passes through air. Wake vortexes from 20 
a lead aircraft 220 may cause the trail aircraft 230 to encoun- 
ter control difficulties and/or turbulence if trail aircraft 230 is 
too close to lead aircraft 220. The size and length of the wake 
vortex created may depend on lead aircraft 220’ s physical 
characteristics (such as size, weight, jet and/or propeller, 25 
wing-span, etc.), and atmospheric conditions (such as visibil- 
ity, wind, temperature, weather, humidity, temperature, pres- 
sure, etc.). 

Similarly, the safe lower limit also depends on the above 
lead aircraft’s 220 physical characteristics and atmospheric 30 
conditions. Additionally, the safe lower limit also depends on 
the trail aircraft’s 230 physical characteristics (e.g., size, 
weight, jet/propeller, wing-span, etc.). 

As an additional embodiment, along the glide path 210, 
there may be a simultaneous runway occupancy (SRO) deci- 35 
sion point where the trail aircraft 230 may execute another 
go-around maneuver 240, 445 (whether it be second attempt, 
third attempt, etc.). Such go-around 240, 445 may be neces- 
sary if a first attempt at the wake vortex decision point is not 
successful. However, when another attempt is made, the go- 40 
around maneuver 240, 445 involves having the trail aircraft 
230 leave the glide path 210. The trail aircraft 230 may then 
circle around and reenter the same or different holding pat- 
tern. Once in a holding pattern, the trail aircraft 230 may 
re-enter the glide path 210 once directed. 45 

As an example, the SRO decision point may be a range of 
points on the glide path 210, such as within one to ten nautical 
miles from the start of runway 260. In another example, the 
SRO decision point may be a specific point on the glide path 
210, such as 5.9 nautical miles from the start of runway 260. 50 
In yet another example, the SRO decision point may be the 
beginning point of any aircraft’s “final approach.” In yet 
another further example, the SRO decision point may be 
measured in time from the runway 260, aircraft’s distance 
between the glide path 21 0 to the runway 260, and/or distance 55 
and/or time between the lead aircraft 220 and the trail aircraft 
230. Also, as another example, the difference in velocity, as 
well as any expected change in velocity, between the lead 
aircraft 220 and the trail aircraft 230 may be considered. 

It should be noted that while the above determinations are 60 
framed in terms of lead aircraft 220 and trail aircraft 230, they 
may also be performed for a series of aircraft, where each trail 
aircraft 230 may be a lead aircraft 220 for the next aircraft to 
enter glide path 210. 

As the above methods and definitions are implementable 65 
by an aircraft landing accident avoidance device 400, refer- 
ence is now made to such device. FIG. 4 shows the various 
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physical modules that make up, and are integrated within, the 
aircraft landing accident avoidance device 400. These mod- 
ules include, but are not limited to, a probability distribution 
function module 405, a safe lower limit module 410, a maxi- 
mum sustainable attempt-to -land-rate module 415, a glide 
slope decision module 420, an actual separation module 430, 
and a warning module 440. 

The probability distribution function module 405 may be 
configured for determining at least two probability distribu- 
tion functions, where one of the probability distribution func- 
tions is runway occupancy time, and where another is landing 
time interval; inter-arrival distance; or a combination thereof; 

The safe lower limit module 410 may be configured for 
determining a safe lower limit on a separation between a lead 
aircraft 220 and a trail aircraft 230 on a glide slope to a runway 
260. The safe lower limit may depend on a multitude of 
factors, such as atmospheric conditions, the physical charac- 
teristics of the lead aircraft 220, and the physical character- 
istics of the trail aircraft 230. Atmospheric conditions may 
include visibility, wind, temp, weather, etc. Physical charac- 
teristics for either the lead aircraft 220 or the trail aircraft 230 
include size, weight, jet/propeller, wing-span, etc. 

The maximum sustainable attempt-to-land-rate module 
415 may be configured for determining a maximum sustain- 
able attempt-to-land-rate (aka target) on the runway 260 
based on the safe lower limit. Like above this rate can be per 
quarter hour, per half hour, per hour, etc. 

The maximum sustainable attempt-to-land-rate module 
415 may also be configured for determining a target separa- 
tion 425. 

The maximum sustainable attempt-to-land-rate module 
415 includes: a maximum sustainable safe throughput mod- 
ule; a maximum sustainable net economic output module; or 
a combination of the above. Sustainable means average. The 
maximum sustainable safe throughput module is configured 
for using a probability of executing the go-around maneuver 
240, 250, 445 and the maximum sustainable attempt-to-land- 
rate. The maximum sustainable net economic output module 
is configured for using a ratio premised on the cost of go- 
around maneuvers 240, 250, 445 to the benefit of landings. 

The glide slope decision module 420 may be configured for 
directing the trail aircraft 230 to enter the glide slope with a 
target separation 425 from the lead aircraft 220 correspond- 
ing to the maximum sustainable attempt-to-land-rate. 

The actual separation module 430 may be configured for 
determining an actual separation 436 between the lead air- 
craft 220 and the trail aircraft 230 while the trail aircraft 230 
is in the glide slope. 

In certain embodiments, a controller device, which may be 
an aircraft instrument, may take and use these stored instruc- 
tions stored in the computer readable storage medium 500 to 
calculate the target separation 425 and actual separation 436. 
In other embodiments, a control tower 270 may use these 
stored instructions to calculate the target separation 425 and 
the actual separation 436. Both the trail separation and the 
actual separation 436 may be designated in distance (e.g., 
nautical miles), time (e.g., seconds, etc.), velocity of each 
aircraft, etc. Where additional multiple instruments (e.g., 
radar, etc.) can determine the target separation and/or actual 
separation distance, such distance may be used to cross check 
the target separation and/or actual separation distance as 
determined via the maximum sustainable attempt-to-land- 
rate module 415 and the actual separation module 430. 

The warning module 440 may be configured for directing 
the trail aircraft 230 to execute a go-around maneuver 240, 
250, 445 if the actual separation 436 approaches the safe 
lower limit at any point in time of an attempted landing. As 
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defined above, “approaches” means the time it takes to reach 
the safe lower limit by taking into consideration the pilot and 
plane’s reaction time — e.g., 1 sec, 5 secs, etc.). The 
go-around maneuver 240, 250, 445 begins at a wake vortex 
decision point. 

Additionally, the warning module 440 considers the reac- 
tion time of the trail aircraft 230 as to when to direct the trail 
aircraft 230 to execute a go-around maneuver 240, 250, 445. 
Reaction time includes, for example, pilot’s reaction time, 
computer response time, mechanical response time, commu- 
nications, sensors, etc. 

It should be noted that while the above determinations are 
framed in terms of lead aircraft 220 and trail aircraft 230, they 
may also be performed for a series of aircraft, where each trail 
aircraft 230 may be a lead aircraft 220 for the next aircraft to 
enter glide path 210. 

As the above methods (and the definitions which the above 
methods rely upon) can also be stored in a computer readable 
storage medium 500, reference is now made to such storage 
medium. Referring to FIG. 5, the computer readable storage 
medium 500 may store instructions for performing a method, 
when executed by a processor or device (such as the aircraft 
landing accident avoidance device 400), for landing aircraft 
on a runway 260. The instructions for such method include: 
determining at least two probability distribution functions 
502, where one of the probability distribution functions is 
runway occupancy time, and where another is: landing time 
interval; inter-arrival distance; or a combination thereof; 
determining a safe lower limit on a separation between a lead 
aircraft 220 and a trail aircraft 230 on a glide slope to the 
runway 260, 505; determining a maximum sustainable 
attempt-to -land-rate on the runway 260 based on the safe 
lower limit 510; directing the trail aircraft 230 to enter the 
glide slope with a target separation 425 from the lead aircraft 
220 corresponding to the maximum sustainable attempt-to- 
land-rate, 515; while the trail aircraft 230 is in the glide slope, 
determining an actual separation 436 between the lead air- 
craft 220 and the trail aircraft 230, 520; and directing the trail 
aircraft 230 to execute a go-around maneuver 240, 250, 445 if 
the actual separation 436 approaches the safe lower limit, 
525, 530. 

Examples of the computer readable storage medium 500 
include, but are not limited to, a compact disc (cd), digital 
versatile disc (dvd), blu ray disc, usb flash drive, floppy disk, 
random access memory (RAM), read-only memory (ROM), 
erasable programmable read-only memory (EPROM), opti- 
cal fiber, electronic notepad or notebook, etc. It should be 
noted that the computer readable storage medium 500 may 
even be paper or other suitable medium in which the instruc- 
tions can be electronically captured, such as optical scanning. 
Where optical scanning occurs, the instructions may be com- 
piled, interpreted, or otherwise processed in a suitable man- 
ner, if necessary, and then stored in computer memory. 

The instructions may be written using any computer lan- 
guage or format. Nonlimiting examples of computer lan- 
guages include Ada, Ajax, Basic, C, C++, Cobol, Fortran, 
Java, Python, XML, etc. 

The instructions on the computer readable storage medium 
500 may be executed by a specific aircraft controller or some 
controller at a control tower 270 that is geared to making safe 
landing determinations and safe (and/or regulated) trailing 
distances between lead aircraft 220 and trail aircraft 230. 

The safe lower limit may depend on a multitude of factors, 
such as atmospheric conditions, the physical characteristics 
of the lead aircraft 220, and the physical characteristics of the 
trail aircraft 230. Atmospheric conditions may include vis- 
ibility, wind, temp, weather, etc. Physical characteristics for 
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either the lead aircraft 220 or the trail aircraft 230 include size, 
weight, jet/propeller, wing-span, etc. 

Once the lower limit has been determined 505, a maximum 
sustainable attempt-to -land-rate can then be determined 510. 
5 Based on the safe lower limit, the maximum sustainable 
attempt-to -land-rate is the target rate at which multiple air- 
craft are allowed to safely land on the same pathway (i.e., 
runway, etc.). The units for such target are time-based, such as 
be per quarter hour, per half hour, per hour, etc. 
to The maximum sustainable attempt-to-land-rate includes 
determining a maximum sustainable safe throughput, a maxi- 
mum sustainable safe net economic output, or a combination 
of the two. Sustainable means average. The maximum sus- 
tainable safe throughput is determined using a probability of 
15 executing the go-around maneuver 240, 250, 445 and the 
maximum sustainable attempt-to-land-rate. The maximum 
sustainable net economic output (e.g., optimality, profit, ben- 
efit, gain, surplus, etc.) is determined using a ratio premised 
on the cost of go-around maneuvers 240, 250, 445 to the 
20 benefit of landings. 

The computer readable storage medium 500 further 
includes instructions for directing the trail aircraft 230 to 
enter the glide slope according to the target separation 425 
after the safe lower limit and the maximum sustainable 
25 attempt-to-land-rate are determined 515. 

After providing such instructions, the computer readable 
storage medium 500 can provide instructions for determining 
the actual separation 43 6 between the lead aircraft 220 and the 
trail aircraft 230 while the trail aircraft 230 is in the glide 
30 slope 520. 

In certain embodiments, a controller device, which may be 
part of the instruments of any aircraft, may take and use these 
stored instructions stored in the computer readable storage 
medium 500 to calculate the target separation 425 and actual 
35 separation 436. In other embodiments, a control tower 270 
may use these stored instructions to calculate the target sepa- 
ration 425 and the actual separation 436. Both the trail sepa- 
ration and the actual separation 436 may be designated in 
distance (e.g., nautical miles), time (e.g., seconds, etc.), 
40 velocity of each aircraft, etc. Where additional multiple 
instruments (e.g., radar, etc.) can determine the target sepa- 
ration and/or actual separation distance, such distance may be 
used to cross check the target separation and/or actual sepa- 
ration distance as determined via the instructions housed in 
45 the computer readable storage medium 500. 

The computer readable storage medium 500 may further 
house instructions for directing the trail aircraft 230 to 
execute a go-around maneuver 240, 250, 445 if the actual 
separation approaches the safe lower limit 525, 530. As 
50 defined above, “approaches” means the time it takes to reach 
the safe lower limit by taking into consideration the pilot and 
plane’s reaction time — e.g., 1 sec, 5 secs, etc.). The 
go-around maneuver 240, 250, 445 begins at a wake vortex 
decision point. 

55 The computer readable storage medium 500 may also fur- 
ther house instructions for directing the trail aircraft 230 to 
attempt a landing if the actual separation does not approach 
the safe lower limit 535. 

Additionally, the computer readable storage medium 500 
60 includes instructions for considering the reaction time of the 
trail aircraft 230 as to when to direct the trail aircraft 230 to 
execute a go-around maneuver 240, 250, 445. Reaction time 
includes, for example, pilot’s reaction time, computer 
response time, mechanical response time, communications, 
65 sensors, etc. 

It should be noted that while the above determinations are 
framed in terms of lead aircraft 220 and trail aircraft 230, they 
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may also be performed for a series of aircraft, where each trail 
aircraft 230 may be a lead aircraft 220 for the next aircraft to 
enter glide path 210. 

The present invention can be applied to any instrumental 
landing and be independent of meteorological conditions. 

Embodiments of Procedures 

I. Determining Probability Distribution Functions 

To determine PDFs using ROT and LTI, IAD, or LTI and 
IAD, mulitlateration data (e.g, from an airport, like Detroit 
Metropolitan Wayne Count airport (DTW)) need to be 
extracted. Out of a possible eighteen fields from multilatera- 
tion data, the present invention may consider five: aircraft 
mode-s, time (t in seconds), longitude (X in meters), latitude 
(Y in meters), and mode-c. The mode-s field is a number of an 
attached transponder that uniquely identifies an aircraft. The 
transponder is generally attached somewhere close to the 
center of the aircraft. The mode-c field is a barometer-based 
value that can be converted to altitude (in feet) by multiplying 
it by 25 and adding 10,000 to the result. However, the 
obtained value is not very reliable for this purpose due to 
pressure change and barometer errors under for different 
weather conditions. Time and position of aircraft are recorded 
every second. 

A. Data Preparation 

The database that stores aircraft landing data at an airport 
may be in Oracle format. SQL+ may be used SQL+ to obtain 
queries. Necessary manipulations and sample extractions 
may be done in MATLAB. 

To start, the following may be sorted in this order: data by 
mode-s and then by time. The time stamp may be changed to 
the format “dd/mm/yyhh:mi:ss .” Basic queries demonstrated 
that the mode-s may be missing for some records. In some 
cases, the mode-s of an entire aircraft track may be missing. In 
other cases, the mode-s of only a few points along a track may 
be missing. Such data points may be eliminated. In the latter 
case, the basic track path may be retained, since the path of the 
aircraft from the other points with mode-s can be linearly 
interpolated. In the former case, the entire track may be dis- 
carded. This action may result in some inter-arrival times that 
are too long. However, because of the available data, losing 
some possible landing records does not have a significant 
impact. Using DTW as the exemplified embodiment, the 
origin (X=0,Y=0) of the Euclidian coordinate may be the 
FAA control tower located between runways 21R and 22L, 
and the Y axis indicates the true north. DTW’s runway 21L, 
and all other runways parallel to it, have a Magnetic angle of 
214.8°. True North and Magnetic North have an angle of 6. 1 ° 
W. Thus, the true angle of runway 21L is 214.8-6.1=208.7°, 
or equivalently 61.3° from the X-axis. Since data may be 
collected in the true coordinates, the same results may be 
ob served by tracking the aircraft course on the runways . In the 
same manner, the true angle of runways 27L/09R and 27R/ 
09L may be calculated as 1.3° from the X-axis. To simplify 
working with the database, the coordinates may be rotated to 
make the runways parallel to the X-axis. To find the aircraft 
position in the rotated coordinates, the observed (XY) posi- 
tion may be multiplied by the rotation matrix R as 


d I Cos (a) Sin(a) (1) 

R l-Sin(a) Cos(a) J 

where a is the rotation angle which is 61.3° for runways 
21L/03R, 21R/03L, 22L/04Rand22R/04L, and 1.3° for run- 
ways 27L/09R and 27R/09L as described before. That is, the 
aircraft position in the rotated Euclidian coordinates is 
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Using the rotation formula (2), the runway coordinates 
may be transformed to the new coordinates. 

Since landings are the subject of study, it is sufficient to 
consider data in a rectangle, which may be referred to as a 
10 “query box.” The sides of the query box are parallel to the 
sides of the runway rectangle, and the box includes the run- 
way and the common landing path extended about 10 nm 
from the runway threshold. For DTW’s runway 21L, for 
example, the rectangle -1350 m<Xr<18500 m, and -3000 
15 m<Yr<-800 m may be considered. Data beyond this rect- 
angle may be dropped from the query. Queries may be 
obtained for every runway for a certain period of time (e.g., a 
week) for data sampling. The time stamp of each of these 
outputs may be transformed to a second-format with respect 
to a time reference. Time zero may be considered as, for 
instance, January 1 . 

Position is recorded at a second rate; however, there may be 
time gaps when position is not recorded. For such cases, if the 
25 gap is at most 1 0 seconds, the timeposition of the aircraft may 
be linearly interpolated between two boundaries of the time 
gap for every second. This interpolation is not applied for 
time gaps greater than 10 seconds. This procedure may be 
implemented in MATLAB. 

30 Other factors that need to be considered include wake 
vortex weight classes, and weather conditions (Instrumental 
Meteorological Condition (IMC) or Visual Meteorological 
Condition (VMC)) information to data records. These may be 
obtained, for example, by matching and search of tables of the 
35 Federal Aircraft Aviation (FAA) aircraft registration data- 
base, including MASTER, ACFTREF, and Aircraft Informa- 
tion tables. The weather condition for every quarter hour may 
be reported in Aviation System Performance Metrics (ASPM) 
database in local time. Considering the time column of the 
40 data, a new column to records may be added to indicate IMC 
and VMC weather condition. 

B. Extracting Samples 

Landings should first be recognized. Then necessary 
records may be extracted through the following steps. 

45 First, for each mode-s, divide all records of a single aircraft 

into separate operations (landings, departures, etc). It is 
assumed that a new operation begins whenever there is a time 
gap of more than 1 5 minutes between any two records of that 
aircraft. Any of these operations may be a landing, departure, 
50 fly over, or a ground operation. 

Second, check if a given operation is a landing on a given 
runway, using DTW’s runways as examples, by checking if it 
passes the following tests: 

Let t min and t max be the first and last times for which the 
55 aircraft is in the “query box.” 

If X(t mi .J-X(t max )>5,000 m, then the aircraft proceeds 
from right to left, and has been long enough in the 
runway direction to be a candidate for a landing on 
runways 21L, 21R, 22L, 22R, 27L, or 27R. Similarly, if 
60 X(t mzw )-X(t max )<-5,000 m, then it is a candidate for a 

landing on runways 03R, 03L, 04R, 04L, 09R, or 09L. 

Check if the aircraft ever crosses the threshold of the spe- 
cific runway and is observed over the runway. 

Third, repeat the second step for all operations and aircraft, 
65 and record their threshold time and location. Record the time 
and location of aircraft when it is first observed outside of the 
runway rectangle after landing, i.e. taxi -in time and location. 
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If the aircraft track disappears over the runway, then exit from 
runway is not recorded, record zero or blank for the exit time. 

Fourth, sort landings in ascending manner, to recognize 
lead aircraft and trail aircraft. Record the location of any trail 
aircraft at the moment its lead aircraft crosses the runway 
threshold. 

Fifth, calculate ROT for any aircraft, and LTI, and IAD for 
any pair of lead aircraft-trail aircraft. 

C. Landing Statistics 

In analyzing system operations, it is important to know the 
peak period for landings. As one embodiment, peak period 
may be defined for a given runway to be a quarter hour with at 
least seven landings on that runway. It may also be important 
to know the proportion of lead aircraft-trail aircraft pairs. The 
aircraft of each pair may be classified into one of these exem- 
plified classes: small, medium, large, B757, or heavy. Any 
combination of two of these classes (e.g., small-large; 
medium-medium; heavy-B757, etc.) may be considered. 

Samples of IAD and samples of LTI may be gathered for 
landings at an airport. Independence of samples may be 
examined by “one-lag scatter plot”. In an embodiment, the 
plot does not demonstrate a specific pattern of dependency 
among the samples; one-lag correlation coefficient is 0.25. 
Higher degrees of lags tend to have lower correlation coeffi- 
cients. Thus, independence of IAD samples, which is 
required for distribution fitting purposes, is accepted. In the 
same manner, independence of LTI samples can be concluded 
by examining the one-lag scatter plot with a related correla- 
tion coefficient of 0.25. 

Based on the samples, histograms and probability distribu- 
tion function (PDF) can be fitted fits for IAD and LTI. For 
practicality, in fitting a distribution, IAD may be limited to a 
minimum of 1.5 nmi. Its distribution may be estimated by 
Erlang where the values represent location (shift), scale, and 
shape parameters, respectively. In one example, where 
samples were based on landings at DTW, Erlang was found to 
be fitted at Erlang(l .5;0.35,6). The mean of the Erlang distri- 
bution is [(location par.)+(shape par.)*(scale par.)], and the 
variance is [(shape par.)*(scale par.)2]. The Maximum Like- 
lihood Estimation (MLE) method may be used for this esti- 
mation and for estimations of LTI and ROT probability dis- 
tributions. The fit may be tested against the Kolmogorov - 
Smirnov test (KS-test) for significance levels less than 0.10. 
The Log-Logistic distribution may provide a slightly better 
fit. In one example, based on the DTW samples, the Log- 
Logistic distribution was found to be Log-Logistic(1.5;1.9, 

4.5) , where values represent location, scale, and shape param- 
eters, respectively. 

When a minimum of 40 seconds is enforced on the DTW 
samples, the probability distribution of LTI may be estimated 
by. Similar to the IAD case, the Log-Logistic (40; 6 1,4. 4) dis- 
tribution may provide a slightly better fit than the Erlang 
distribution, which is a specific case of the gamma distribu- 
tion. The Erlang fit is accepted by KS-test for significance 
levels of 0.05 or smaller. 

Samples (e.g., 669 samples at DTW) of ROT in peak IMC 
periods may be obtained. Those collected at the DTW 
example may be concluded to be independent because the 
N-lag correlation coefficient, N=1 , 2, . . . , is less than or equal 
to 0.08, and because the one-lag scatter plot does not show 
any specific pattern of relationship. The distribution of ROT 
may be estimated using three different distributions — 
gamma, beta, and normal. Using the MLE method, the best 
fits for these samples were determined to be Gamma(25;2.8, 

8.5) in the enforced range of (25, oo) s, Beta(25,110;6.1,15.4) 
in the enforced range of (25, 1 1 0) s, where 3rd and 4th values 
represent shape parameters, andN(49,8.12) in the open range 


of (- 00 , 00 ). Gamma is the best fit among these three with the 
maximum likelihood criterion. However, the beta distribution 
might be preferred because, as in real situations for ROT, it 
has lower and upper bounds (for example, it can not be nega- 
5 tive). The normal distribution is rejected for ROT samples in 
the 0.1 significance level. Mean and variance of Beta(L,U:a, 
(3) are: 


to 



(a + f3) 2 (a + f3 + 1 )‘ 


(3) 
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As another embodiment, it would be helpful to know if 
ROT is different under IMC and VMC weather conditions. 
Also of interest is the probability (or frequency) that the LTI 
between two consecutive aircraft is less than the ROT of the 
20 leading aircraft. This probability may be represented by 

P{LTI Ar ^ +1 <ROT Ar }, k=l, 2, It may be named as “runway- 

related approach risk.” 

Empirical and a theoretical “point estimation” for 
P { YI\ kk+ 1 <ROT jt } , k=l, 2, . . . , in peak periods may be 
25 conducted for pairs of aircraft with a certain separation stan- 
dard (e.g., 3 nmi). 

Empirical Method Where a 45 degree line, for example, is 
plotted on an LTI (x-axis) and ROT (y-axis) grid, and points 
lie above such line, these points represent events where 
30 LTI ArA . +1 <ROT Ar . In one study involving data collected (e.g., 
4,312 landings) at DTW, the frequency of LTI^ +1 <ROT A , for 
k=l, . . . , 4312, is 0.0007 — that is, 3 out of 4313 landings, 
assuming that no such event occurred in non-peak periods. 

Theoretical Method Using the sample data collected (e.g., 
35 4,3 1 2 landings) at DTW, the probability distribution fits was 
calculated as ROT~Beta(6.1,14.5) in the range (25,1 10), and 
LTI-Erlang(40;l 1,6) to estimate P{LTI<ROT}. FIG. 7 shows 
an overlap of the probability of the distributions of ROT and 
LTI. Because there is an overlap between LTI and ROT, then 
40 P{LTI<ROT} is positive. In fitting the PDF for LTI, samples 
of LTI that were not considered include those in which the 
corresponding ROT could not be obtained. Let g ROI (*) rep- 
resent PDF of ROT, and F i77 (*) represent Cumulative Density 
Function (CDF) of LTI. Then, 

45 


50 


p = {LTI < ROT} 



< ROT | ROT = x}-g R0T (x)dx 



p{LTI < x}-g R0T (x) dx 


F LTl {x)-g R0 T{x)dx. 


(4) 


55 

II. Risk Free Landing (Enforcing a Go-around Procedure) 
Generally, it is desired that the chance of an SRO or a WV 
hazard, whether moderate or severe, be nearly or exactly zero. 
Conventionally, increasing the target separation between suc- 
60 cessive aircraft decreases these two risks. The risk can also be 
reduced by implementing go-around procedures. For 
example, if two aircrafts will be on the runway at the same 
time, the trailing aircraft can execute a go-around procedure 
to avoid a SRO. In reality, the go-around is not always taken. 
65 In these exemplified embodiments, it may be assumed that 
an aircraft is always enforced to execute a go-around when- 
ever separation minima are not or will not be met. In addition, 
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perfect information may be assumed in making this decision. 
With these assumptions, the risk of an SRO or a wake-vortex 
hazard is exactly zero, though there is possibly an increase in 
the number of go-arounds. 

Making the system risk-free by enforced GA creates a 5 
different dynamic and may change the optimal level of opera- 
tions (i.e., the best number of attempts per hour). The optimal 
level of attempts per hour depends on the GA probability 
P{GA} . This section calculates this probability for two cases 10 
of with and without wake vortex effect. 

In the approach/landing process, two different aircraft 
flows can be recognized: the flow through the glide slope go 
measured in landing attempts per hour, and the flow through 
the runway (or simply throughput) X measured in actual land- 1 5 
ings per hour. FIG. 6 demonstrates this dynamic with 
enforced GA procedures. When the following aircraft is at 
decision point 2 (DP2) (e.g., 8 nmi from runway threshold), 
the controller/pilot decide(s) whether or not to take a GA 
procedure to avoid the risk of encountering a hazardous wake 20 
from the leading aircraft. If the separation is less than a 
specific value, x 0 , at this point, then the following aircraft 
must go -around (GA) to a holding position and return to the 
glide slope when cleared to attempt again. Such a minimum 
WV safe separation exists and can be estimated using wake 25 
vortex theories, and/or field observations. This operation may 
be referred to as the “wake vortex GA” or “wake vortex 
missed approach” procedure in contrast with the well known 
GA procedure executed to avoid a SRO. The latter may be 
referred to as a “SRO GA” or “SRO missed approach”. For 30 
illustrative purposes, DP2 is considered to be 8 nmi from the 
threshold, and x 0 =65 s separation as the minimum WV safe 
separation. If a safe separation is achieved at DP2, then the 
aircraft can continue the approach. At DPI, which is called 
the decision height, the follower decides whether or not to 35 
execute a go -around to avoid simultaneous runway occu- 
pancy with the leading aircraft. 

T the total GA probability can be defined as p. It should be 
noted that p is a function of go, the number of attempts per 
hour. The average GA rate (number of go-arounds per hour) is 40 
P(go)*go and the average successful landing rate is X(Go)=[l-p 
(go)] - go . The rate of aircraft attempting to land is the arrival rate 
of aircraft Z(go) plus the rate of aircraft executing a go-around 
P(go) go. As a check for consistency, the attempt rate is [1-p 
(go)]-go+p(go)-go which equals go. 45 

In addition, the following assumptions may be made: 

YI\ k+l k and ROT^ are independent random variables 

The separation is minimized at DP2 and remains 
unchanged afterwards until the touchdown. In other 50 
words, the separation at DP2 equals LTI 

Shifting LTI to the right or left does not change the shape of 
its probability distribution 

Zero risk assumed for execution of both GA procedures 55 

GA are absolutely respected and enforced at both decision 
points 

The number of GA in an hour is not restricted 

Wake vortex GA and SRO GA conditions never simulta- 60 
neously occur for a pair. That is, no simultaneous go- 
around for aircrafts k and k+1 happens for all 
k=l, 2, . . . . 

where LTI^ +1 is the landing time interval between aircraft k 
and k+1 measured at the runway threshold (in seconds), and 65 
ROT k is the runway occupancy time of aircraft k (measured in 
seconds). 
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A. Go-Around Probability Assuming No Wake Vortex 
Effect 

In this section, the possible go-around at DP2 is ignored. In 
other words, the only tiling considered is the risk of a SRO. 
The risk of a hazardous wake vortex encounter is not consid- 
ered. The probability of a SRO is 


P{SRO} = P{LTI < ROT & Follow aircraft lands} 

= P{Follow aircraft lands | LTI < ROT} * 

P{LTI < ROT}. 

This probability can be reduced to zero by enforcing the 
go-around procedure. In this case, 

Pj Follow aircraft landslLTI<ROT}=0, and 

Pi^^jGAj^lLTKROT}. (5) 

Probability distribution functions of peak period LTI and 
ROT are estimated for DTW for the pairs of FAA 3 nmi 
minimum separation pairs. These follow-lead aircraft pairs 
include S-S, L-S, B757-S, H-S, L-L, B757-L, and H-L. The 
estimations are used here for methodology illustrations. LTI 
is the peak period distribution calculated for arrival of aircraft 
to the glide slope (or the final approach fix) with the rate go. 

FIG. 7 illustrates the LTI and ROT probability distributions 
obtained for 3 nmi pairs. p(go) is estimated as 0.004, the mean 
of LTI is 106 s, and the average number of attempts per hour 
(during peak periods) is go=3 600/ 106=34 attempts/h. In this 
period, no go-around was observed, so go=Z, but 
P{SRO}=0.004 instead of 0.0. 

Changes in arrival rates (that is, changes in A, or go) can be 
modeled by shifting the LTI distribution to the left or right 
(i.e., by changing the location parameter). Shifting the LTI 
distribution also changes the probability that LTI<ROT, or 
equivalently, the probability of a go-around. 
P(go)=P{LTI<ROT } as indicated by the broken line in FIG. 8. 
Shifting LTI toward left or right (i.e., changing the location 
parameter) provides different values for go and p (go) assuming 
other parameters of the distribution remain unchanged by 
shifting. p(go)=P{LTI<ROT} is given by broken line in FIG. 
8 . 

B. Total go-Around Probability with Wake Vortex Effect 

When the wake vortex effect is taken into account and WV 
GA procedure is in place, an aircraft would possibly miss the 
approach for two reasons at two different points: WV safe 
threshold and runway safe threshold (decision height). For 
this situation, let x 0 to be the minimum wake vortex safe 
separation of successive aircraft given in seconds, and L to be 
the lower limit of LTI distribution; e.g., L=40 s. To calculate 
P(GA) two cases for L and x 0 can be considered. 

Case 1 : L<x0. For this case, 


P(GA)= F LTI (y)dF R oT(y) + F LTI (xo)F R oT(xo), 

JXQ 

where F ir/ and F^ or are CDF of LTI and ROT, respectively. 
Detailed calculations may be given as follows: For the case 
L<x 0 , note that P(GA)=P{ separation at DP2<x 0 or, exclu- 
sively, separation at DPl<ROT}=P{LTI<x 0 or LTI<ROT}. 
Then, 

J P(GA)= J PLTI<JC 0 )+ J P{LTI<ROT and LTI^x 0 }. (1. 1) 
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On the other hand, 


P{LTI < ROT and LTI>x 0 }= f P 

JXQ JXQ 
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( 1 . 2 ) 
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FlTI (*0 ) + F LT i {Xq)F R qt (Xq ), 


where F LTI and F ROT are CDF of LTI and ROT, respectively. 
Joint distribution of LTI and ROT is broken into multiplica- 25 
tion of their marginal distributions because of their indepen- 
dence. Plugging in (1.2) in (1.1) gives equation (6). 

Case 2: L>=x 0 . This case means that LTI shifted to the right 
as much that its lower point L is above the wake vortex safety 
threshold of x 0 . No wake vortex GA would ever occur in this 30 
case, and P{GA} is obtained from equation (5). 

Total p(oo)=P{GA} is shown in FIG. 8 for peak period 
DTW IMC distributions, which are given in FIG. 7, and for 
x 0 =65 s. This figure also shows p l and derivative of p(oo) in 
terms of oo. The derivative dp/doo, which will be used in the 35 
following sections, is multiplied by 1 0 to make it more visible 
in this figure. 

III. Maximizing Runway Throughput 
The present invention looks to find the relation between 
glide-slope rate of to attempts/h, runway throughput X(co), 40 
and go around probability p(oo). Then one can find out for 
what values of oo and p throughput X is maximum, while 
P{ SRO}=0 is maintained by enforced GA procedure. In other 
words, the objective is maximization of the runway through- 
put. That is, 45 

Maximize X(co)=[l-p((o)]-co. (7) 

This model is the same for both with and without WV effect 
assumptions. However, p(oo) differs depending on each of 
these cases, as discussed below. 50 

A. Maximum Throughput without Wake Vortex Effect 
In this case p(oo) is calculated from equation (5) and 
plugged into the problem (4) to maximize the throughput. 
FIG. 9 provides L(oj) in the left axis for the distributions in 
hand. X(oo) is calculated for all pairs of (co, p(co)) in the right 55 
axis. By increasing the rate of attempts co, the percentage of 
GA increases but the percentage of successful landings 
decreases. After a point, the decrease in the rate of successful 
landings dominates the increase in the rate of attempts. In 
other words, throughput L(oo) has a unique maximum or 60 
optimal point. This scenario can also be explained in math- 
ematical terms. p(oo)>0 is increasing, and l-p(oo) is decreas- 
ing in co. So after a point, decrease of l-p(co) dominates 
increase of co and [l-p(co)]*co would have a maximum. 

For distributions in hand, the optimal (co,X,p) is (46.5, 39.6, 65 
0.148), as indicated in FIG. 9. To have a stable system, the 
arrival rate to the TRACON, X, is adjusted so that co is main- 
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tained in the optimal level of 46.5 attempts/h. 39.6 landings/h 
is the maximum and optimal throughput. 

B. Maximum Throughput with Wake Vortex Effect 

The previous section analyzed the optimal level of landing 
attempts when the wake vortex constraint was relaxed and 
SRO was the only risk factor. In this section, wake vortex risk 
is also considered in maximization of runway throughput. 
Using P{GA}=p(oo) from equations (6) and (5) in problem 
(3), the following solution may be obtained. 

Throughput X(co) in terms of co attempts/h is shown in FIG. 
10 with the solid line (indicated by r=0). Other plots of this 
figure are the subject of next section. The optimal (co,X,p) is 
about (37.0, 33.5, 0.08) and co-p(co) is 3 go-around per h. 

WV is costing the system 39.6-33.5=6 landings per peak 
hour. Thus WV is costing the system about 21,900 landings 
per year if there is 10 hours of peak period every day. This 
scenario is for the cases of 3 nmi pairs, in which the majority 
of them are large aircrafts. Now this value can be multiplied 
by R to obtain dollar value WV cost estimation. For R=$ 1 000, 
it is about $22 m per year for a single landing runway under 
IMC. 

IV. Maximizing Net Economic Gain 

Maximizing the number of (successful) landings does not 
necessarily guarantee the overall economic optimality of the 
landing operations. The reason is because costs and gains 
from the operations are important optimality parameters and 
will be taken into account. For the landing operations when 
absolute safety is guaranteed by enforced GA procedure, the 
economic profit/benefit to the regional/whole economy (in- 
cluding airlines, airport, employees, etc) is the result of a 
successful landing. The overall cost to the regional/whole 
economy associates with the go-around procedure. The costs 
of successful landings are embedded in the landing profits, 
that is, total revenue to all stakeholders minus operational 
cost, except the cost of GA execution. The net economic gain 
or surplus, that is, total gain minus total cost of an hour of peak 
period operations, is desired to be maximized with respect to 
the number of attempts/h rate, to. Since this net economic gain 
is a random variable, we consider maximizing its expected 
value, ES=E{economic surplus}. 

The gain from one successful landing is R which occurs 
with probability 1 -p(co) for every landing attempt. The loss of 
one landing attempt is the cost of go-around C which occurs 
with probability p(co). Thus, since the number of attempts per 
hour is on, then the expected value of the net gain from hourly 
landing attempts is ES(co;R,C) given in (4) and the optimiza- 
tion objective is 

Maximize C)=vy[(l-p(<d)yR-p((ayC] (8) 

Some dollar values for R and C can be considered to 
illustrate economic behavior of the system. For any given 
type of aircraft, C is the summation of cost components such 
as passenger delay cost, disturbed schedules cost of down- 
stream flights, take off cost, aircraft operations cost, and 
airport cost. Any of these cost components depend on param- 
eters such as aircraft load factors and the arrival rate at a given 
time, which are uncertain. Thus, C is a random variable. 
However, its expected value is considered as a suitable esti- 
mation of this parameter. Estimation of R and C is part of on 
ongoing research. However, at the time being and for the sake 
of illustration, it is considered that C=$4,000 for a large 
aircraft in a peak hour. Three scenarios of $ 1 ,000, $2,000, and 
$4,000 are considered for R here. 

A. Economic Optimality without Wake Vortex Effect 

For this case p(co) is calculated from equation (5), and fed 
into problem (4). FIG. 11 is a plot of ES(oo;R,C) in thousands 
of dollars, for assumed values of R and a fixed C=$4,000. 
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To obtain a more general form, ES can be written in terms 
of the ratio C/R. Factoring out R in equation (6) gives 

ES(ar, R, C) = R ■ to ■ [(1 - pM) - pM ■ |] 


Thus, ES(oo;R,C) is a multiplication of constant dollar 
value R and a function of co. Define the latter function to be 
g( # ) as follows: 

g-((o;r)=( 0 - [1 -(1 +r)/?((o)] 
or 

g((o;r)=A,(<i ))- r -co -p (to) , 

where r=C/R. Thus, maximizing g(co;r) is equivalent to maxi- 
mizing ES(oo;R,C). So the problem reduces to the more gen- 
eral form of 

Maximize g(co; r) = co ■ [1 - (1 4 - r)p(a >)] (9) 

= A(w) - r- co- p(co). 

FIG. 12 illustrates g(co;r) forr=0, 1, 2, 4 and co in [28,55]. 

Note that the optimal solution co !is (r) r =Argmax{g((o;r)} 
depends on the ratio r=C/R. The derivative of g(oo;r) with 
respect to co is zero at co*. ES(co;R,C) and g(co;r) have the 
following interesting properties (at least for the LTI and ROT 
distributions in hand): 

Property 1 : For given R and C, ES(co;R,C) and g(co;r) are 
unimodal in a practical peak period rates of attempts/h. That 
is, they have a unique maximum in this range. This is the 
necessary optimality condition. This fact can be justified by 
visual investigations of FIGS. 9 and 10 for given DTW peak 
period landing distributions. However, this fact remains to be 
proven mathematically. In FIG. 9 , ES(co;R,C) has a unique 
maximum for any given R and C. In FIG. 10 , g(co;r) has a 
unique maximum for any given r. 

Property 2: g(co;r) decreases as r increases for any fixed co. 
This is seen from equation (9). The only term that includes r 
is the term inside the brackets which is decreasing in r. 

Property 3: g(co;r) ^X(co)=(l-p(co))*co whenr— >0 + , which 
is obvious from equation (9). So, g(co;r) is bounded by 
X(co)=[l-p(co)]*co and is below it, see FIG. 12. 

Property 4: co*(r)=Argmax{g(co;r)} is decreasing in r for 
2 8<co<Argmax{ dp/ dco } . 

The proof for Property 4 is by contradiction. co*(r) is the 
point where dg/dco=0 or p(w*)+co*-p , (co*)=l/(l+r). Increas- 
ing r decreases the right hand side and consequently the left 
hand side of this equation. On the other hand, co, p(co), and 
dp/dco are all increasing in co for co<Argmax{dp/dco}. Thus if 
co does not decrease, the LHS will not decrease. This showing 
is a contradiction, and completes the proof. 

These properties are seen in FIGS. 11 and 12 . They are 
unimodal (property 1). co* decreases as r increases (property 
2). Argmaxj dp/dco} is calculated 49.7 attempts/h for the LTI 
and ROT in hand. g(co;r) is below X(oo). Peak of g(co;r) moves 
down and left by increasing the relative penalty of go -around 
to the landing benefit. Properties of g(co;r) imply that the 
highest value for the optimal number of attempts per hour and 


18 

the upper bound of the optimal throughput is the maximal 
point of X(co). This scenario is achieved when C is much 
smaller than R. 

5 In such a case, the problem reduces to maximizing X(co)= 
[l-p(co)]*co scenario, as discussed above. The highest 
throughput value is a good estimation of the average runway 
landing capacity. In the example provided, this capacity is 
39.6 landings/h for the case that wake vortex effect is ignored, 
10 and GA procedure is enforced. 

As examples, forr=l, (co, X, p)*=(40.0, 38.2, 0.045). For 
r=2, (co, X, p)*=(38, 37.1, 0.024). The latter means that to 
maximize the expected value of the net economic gain (sur- 
plus) from the landing operations, when go-around cost C is 
2 times larger than landing profit R, the average glide slope 
throughput shall be adjusted at 38.2 attempts/h which gives 
37.0 successful landings/h. Note that in this case, system 
throughput is 3 landings/h more than the current level of 34 
20 landing/h with associated P{SRO}=0.004, as indicated in 
FIG. 9. 

B. Economic Optimality with Wake Vortex Effect 

In this case the problem is maximizing (6), or equivalently 
25 (4), where P{GA}=p(co) is calculated from equations (6) and 
(5). Note that oo=3600/mean(LTI) where LTI is in seconds. 
All properties of g(co;r) in the section above are still valid for 
this new p(co). Justification procedures are the same that are 
provided above. For DTW peak period IMC distributions of 3 
30 nmi pairs derivative of p(co) is maximized at 41.3 attempts/h 
when wake vortex safe threshold is 65 s or 2.2 nmi, as seen in 
FIG. 8. This is the condition for property 4. 

The optimal results for DTW distributions are provided in 
TABLE 1 and FIG. 10 for x 0 =65, 70, and 75 seconds. FIG. 13 
35 is visualization of TABLE 1 . The optimal solution (co, X, p)* 
is (36.8, 33.6, 0.087) for (x 0 ,r)=(65, 0). Since r=0, that is, cost 
of go around is negligible in comparison with the landing 
profit, then this is the optimal number of landings/h. So the 
4Q average landing capacity of the system is 33.6 landings/h 
independent of the market condition, (co, X, p)* is (32.7, 32.3, 
0.014) for (x 0 ,r)=(65,2), meaning that the optimal throughput 
is 32.3 landings/h if C=2R in the market, that is, the cost of 
go-around is two times bigger than the profit gained from a 
45 successful landing. Note that in TABLE 1, optimal through- 
put decreases as safe WV threshold increases. However, 
eventually the safe threshold is a certain number and once it is 
recognized other cases become irrelevant. 

Some capacity estimations would be helpful. In present 
50 literature, the reciprocal of the minimum safe separation is 
sometimes considered as an estimation of the capacity. Using 
that known method for the safe minimum separation of 65 s, 
one obtains the capacity of 55.4 landings per hour. One can 
55 intuitively recognizes that this number is too high since in 
practice the capacity is generally between 30 and 40 landings/ 
h. The problem with this method is that it ignores the proba- 
bilistic nature of the process. Achieving this level of through- 
put requires that the mean of LTI to be adjusted at 65 s which 
60 implies P{LTI<ROT}=0.30 from FIG. 9. In other words, with 
the enforced go-around procedure, there will be more than 
30% loss of attempted landings at DPI . Such case would lead 
to the throughput level of less than 32 landing/h as can be seen 
65 in FIG. 10. Further, in such an operation the system should 
tolerate the high cost of 1 3 go-around/h. So this known recip- 
rocal method is not suitable. 
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TABLE 1 


Optimal Values for Different safe WV threshold and r 

WV 

threshold 

r 

(0* 

X* 

p* 

x 0 = 65 s 

0 

36.8 

33.6 

0.087 


1 

33.7 

32.8 

0.026 


2 

32.7 

32.3 

0.014 


3 

32.1 

31.9 

0.009 


4 

31.8 

31.6 

0.007 

x 0 = 70 s 

0 

34.7 

32.0 

0.079 


1 

32.1 

31.4 

0.024 


2 

31.3 

30.8 

0.013 


3 

30.8 

30.5 

0.009 


4 

30.4 

30.2 

0.006 

x 0 = 75 s 

0 

33.0 

30.5 

0.073 


1 

30.7 

30.0 

0.022 


2 

29.9 

29.5 

0.012 


3 

29.5 

29.2 

0.008 


4 

29.1 

29.0 

0.006 


The present invention’s methodology, which considers the 
probabilistic behavior of the system, estimated the average 
capacity of 33.6 landings per hour for 65 s safe WV threshold 
and risk free landings. Economic considerations may reduce 
the optimal throughput to about 32 landings/h with maxi- 
mized net economic gain, 1% go-around, and risk free (safe) 
operations. 

Additionally, the process may determine a throughput that 
maximizes an economic throughput, based on the cost-ben- 
efit ratio of landings to go-around maneuvers. The benefit of 
a successful landing may include the net benefit to one or 
more of the airlines, passengers, airports, employees, etc. In a 
further embodiment, the net benefit of a successful landing is 
the total revenue to all beneficiaries minus operational costs, 
except the costs of a go-around maneuver. The cost of go- 
around maneuvers may include the cost to one or more of the 
airlines, passengers, airports, employees, etc. in lost time, 
fuel, goodwill, salaries, increased wait time for other aircraft 
in the holding pattern, etc. In one embodiment, the cost of a 
go-around maneuver and the benefit of a successful landing 
are used to help determine economically the most beneficial 
number of attempted landing and projected successful land- 
ings versus go-around maneuvers. In a further embodiment, 
the economically most beneficial number of attempted land- 
ing sets a lower limit on the separation distance between 
aircraft to avoid the risk of SRO and/or WV. 

Many of the elements described in the disclosed embodi- 
ments may be implemented as modules. A module (some- 
times referred to as element, component, or mechanism) is 
defined here as an isolatable element that performs a defined 
function and has a defined interface to other elements. The 
modules described in this disclosure may be implemented in 
hardware, software, firmware, wetware (i.e., hardware with a 
biological element) or a combination thereof, all of which are 
behaviorally equivalent. For example, modules may be 
implemented as a software routine written in a computer 
language (such as C, C++, Fortran, Java, Basic, Matlab or the 
like) or a modeling/simulation program such as Simulink, 
Stateflow, GNU Octave, or Lab VIEW MathScript. Addition- 
ally, it may be possible to implement modules using physical 
hardware that incorporates discrete or programmable analog, 
digital and/or quantum hardware. Examples of program- 
mable hardware include: computers, microcontrollers, 
microprocessors, application-specific integrated circuits 
(ASICs); field programmable gate arrays (FPGAs); and com- 
plex programmable logic devices (CPLDs). Computers, 
microcontrollers and microprocessors are programmed using 
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languages such as assembly, C, C++ or the like. FPGAs, 
ASICs and CPLDs are often programmed using hardware 
description languages (HDL), such as VHSIC hardware 
description language (VHDL) or Verilog, that configure con- 
5 nections between internal hardware modules with lesser 
functionality on a programmable device. Finally, it needs to 
be emphasized that the above mentioned technologies are 
often used in combination to achieve the result of a functional 
module. 

10 While various embodiments have been described above, it 
should be understood that they have been presented by way of 
example, and not limitation. It will be apparent to persons 
skilled in the relevant art(s) that various changes in form and 
detail can be made therein without departing from the spirit 
15 and scope. In fact, after reading the above description, it will 
be apparent to one skilled in the relevant art(s) how to imple- 
ment alternative embodiments. Thus, the present embodi- 
ments should not be limited by any of the above described 
example embodiments. In particular, it should be noted that, 
20 for example purposes, the above explanation has focused on 
the example of a single WV decision point and a single SRO 
decision point. However, one skilled in the art will recognize 
that additional embodiments could include more than one 
WV decision point and/or SRO decision point, or a con- 
25 tinuum of such points on the glide slope. 

In addition, it should be understood that any figures which 
highlight the functionality and advantages, are presented for 
example purposes only. The disclosed architecture is suffi- 
ciently flexible and configurable, such that it may be utilized 
30 in ways other than that shown. For example, the steps listed in 
any flowchart may be re-ordered or only optionally used in 
some embodiments. 

Further, the purpose of the Abstract of the Disclosure is to 
enable the U.S. Patent and Trademark Office and the public 
35 generally, and especially the scientists, engineers and practi- 
tioners in the art who are not familiar with patent or legal 
terms or phraseology, to determine quickly from a cursory 
inspection the nature and essence of the technical disclosure 
of the application. The Abstract of the Disclosure is not 
40 intended to be limiting as to the scope in any way. 

Finally, it is the applicant’s intent that only claims that 
include the express language “means for” or “step for” be 
interpreted under 35 U.S.C. 112, paragraph 6. Claims that do 
not expressly include the phrase “means for” or “step for” are 
45 not to be interpreted under 35 U.S.C. 112, paragraph 6. 

What are claimed are: 

1. A method for landing aircraft on a runway using an 
aircraft landing accident avoidance device, comprising: 

50 a. determining at least two probability distribution func- 
tions, where one of the probability distribution functions 
is runway occupancy time, and where another is: 

1. landing time interval; 

2. inter-arrival distance; or 

55 3. a combination thereof; 

b. determining a safe lower limit on a separation between a 
lead aircraft and a trail aircraft on a glide slope to the 
runway; 

c. determining a maximum sustainable safe attempt-to- 
60 land rate on the runway based on the safe lower limit and 

the probability distribution functions; 

d. directing the trail aircraft to enter the glide slope with a 
target separation from the lead aircraft corresponding to 
the maximum sustainable safe attempt-to-land rate; 

65 e. while the trail aircraft is in the glide slope, determining 
an actual separation between the lead aircraft and the 
trail aircraft; and 
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f. directing the trail aircraft to execute a go-around maneu- 12. The device according to claim 11, wherein the maxi- 

ver if the actual separation approaches the safe lower mum sustainable safe throughput module is configured for 

limit . using a probability of executing the go -around maneuver and 

2. The method according to claim 1, wherein thego-around the maximum sustainable safe attempt-to-land rate. 


maneuver begins at a wake vortex decision point. 

3. The method according to claim 1, wherein determining 
a maximum sustainable safe attempt-to-land rate includes 
determining: 

a. a maximum sustainable safe throughput; 

b. a maximum sustainable safe net economic output; or 

c. a combination of the above. 

4. The method according to claim 3, wherein the maximum 
sustainable safe throughput is determined using a probability 
of executing the go-around maneuver and the maximum sus- 
tainable safe attempt-to-land rate. 

5. The method according to claim 3, wherein the maximum 
sustainable net economic output is determined using a ratio 
premised on the cost of go-around maneuvers to the benefit of 
landings. 

6. The method according to claim 1, wherein the safe lower 
limit depends on atmospheric conditions. 

7. The method according to claim 1, wherein the safe lower 
limit depends on the physical characteristics of the lead air- 
craft and the physical characteristics of the trail aircraft. 

8. The method according to claim 1, wherein the timing of 
when to direct the trail aircraft to execute a go-around maneu- 
ver considers the reaction time of the trail aircraft. 

9. An aircraft landing accident avoidance device compris- 
ing: 

a. a probability distribution function module configured for 
determining at least two probability distribution func- 
tions, where one of the probability distribution functions 
is runway occupancy time, and where another is: 

1 . landing time interval; 

2. inter-arrival distance; or 

3. a combination thereof; 

b. a safe lower limit module configured for determining a 
safe lower limit on a separation between a lead aircraft 
and a trail aircraft on a glide slope to a runway; 

c. a maximum sustainable safe attempt-to-land rate module 
configured for determining a maximum sustainable safe 
attempt-to-land rate on the runway based on the safe 
lower limit and the probability distribution functions; 

d. a glide slope decision module configured for directing 
the trail aircraft to enter the glide slope with a target 
separation from the lead aircraft corresponding to the 
maximum sustainable safe attempt-to-land rate; 

e. an actual separation module configured for determining 
an actual separation between the lead aircraft and the 
trail aircraft while the trail aircraft is in the glide slope; 
and 

f. a warning module configured for directing the trail air- 
craft to execute a go-around maneuver if the actual sepa- 
ration approaches the safe lower limit. 

10. The device according to claim 9, wherein the go-around 
maneuver begins at a wake vortex decision point. 

11 . The device according to claim 9, wherein the maximum 
sustainable safe attempt-to-land rate module includes: 

d. a maximum sustainable safe throughput module; 

e. a maximum sustainable net economic output module; or 

f. a combination of the above. 


5 13. The device according to claim 11, wherein the maxi- 

mum sustainable net economic output module is configured 
for using a ratio premised on the cost of go-around maneuvers 
to the benefit of landings. 

14. The device according to claim 9, wherein the safe lower 
to limit depends on atmospheric conditions. 

15. The device according to claim 9, wherein the safe lower 
limit depends on the physical characteristics of the lead air- 
craft and the physical characteristics of the trail aircraft. 

16. The device according to claim 9, wherein the warning 
15 module considers the reaction time of the trail aircraft as to 

when to direct the trail aircraft to execute a go-around maneu- 
ver. 

17. A computer readable storage medium including 
instructions for performing a method, when executed by a 

20 processor, for landing aircraft on a runway using an accident 
avoidance processor, the method comprising: 

a. determining at least two probability distribution func- 
tions, where one of the probability distribution functions 
is runway occupancy time, and where another is: 

25 1. landing time interval; 

2. inter-arrival distance; or 

3. a combination thereof; 

b. determining a safe lower limit on a separation between a 
lead aircraft and a trail aircraft on a glide slope to the 

30 runway; 

c. determining a maximum sustainable safe attempt-to- 
land rate on the runway based on the safe lower limit and 
the probability distribution functions; 

d. directing the trail aircraft to enter the glide slope with a 

35 target separation from the lead aircraft corresponding to 

the maximum sustainable safe attempt-to-land rate; 

e. while the trail aircraft is in the glide slope, determining 
an actual separation between the lead aircraft and the 
trail aircraft; and 

40 f. directing the trail aircraft to execute a go-around maneu- 
ver if the actual separation approaches the safe lower 
limit. 

18. The computer readable storage medium according to 
claim 17, wherein the go-around maneuver begins at a wake 

45 vortex decision point. 

19. The computer readable storage medium according to 
claim 17, wherein determining a maximum sustainable safe 
attempt-to-land rate includes determining: 

a. a maximum sustainable safe throughput using a prob- 
50 ability of executing the go-around maneuver and the 

maximum sustainable safe attempt-to-land rate; 

b. a maximum sustainable net economic output using a 
ratio premised on the cost of go-around maneuvers to the 
benefit of landings; or 

55 c. a combination of the above. 

20. The computer readable storage medium according to 
claim 17, wherein the safe lower limit depends on: 

a. atmospheric conditions; 

b. the physical characteristics of the lead aircraft; and 
60 c. the physical characteristics of the trail aircraft. 



